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Abstract. A model of compact object coupled to inhomogeneous anisotropic dark 
energy is studied. It is assumed a variable dark energy that suffers a phase transition 
at a critical density. The anisotropic A— Tolman-Oppenheimer-Volkoff equations are 
integrated to know the structure of these objects. The anisotropy is concentrated on 
a thin shell where the phase transition takes place, while the rest of the star remains 
isotropic. The family of solutions obtained depends on the coupling parameter between 
the dark energy and the fermionic matter. The solutions share several features in 
commmon with the gravastar model. There is a critical coupling parameter that gives 
non-singular black hole solutions. The mass-radius relations are studied as well as the 
internal structure of the compact objects. The hydrodynamic stability of the models 
is analyzed using a standard test from the mass-radius relation. For each permisible 
value of the coupling parameter there is a maximum mass, so the existence of black 
holes is unavoidable within this model. 
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1. Introduction: 

A star that consumed its nuclear fuel, if massive enough, will end its life as a compact 
object or a black hole. The black holes described by the Einstein theory of gravity 
contain singularities. However quantum effects must be taken into account at high 
curvature values, or short distances compared with the Planck length scale pQ. The 
quantum effects can render an inner portion of the black hole as a de Sitter spacetime. 
This was considered by several researchers, see for example: Gliner [2], Gliner and 
Dymnikova [3], Poisson and Israel [4J, and references therein. Recently, Nicolini et al. 
[7j, [8], obtained regular black hole solutions with an inner de Sitter portion through a 
non-commutative geometric approach to quantum gravity. 

Analogous concepts were extended to a set of models which replace the whole black 
hole region with a de Sitter spacetime [5], [6]. In fact, Chapline et al [5], built a model 
of dark energy stars based on the analogy between a superfluid condensate near its 
critical point with the neighborhood of an event horizon. The dark energy stars have 
a de Sitter spacetime matched with a Schwarzschild exterior spacetime. The surface of 
phase transition is closely located above the Schwarzschild radius. The model of Mazur 
and Mottola (MM) [6] is conceptually different from the dark energy stars, but possess 
several common features inspired in it. Their model is called "gravitational vacuum 
stars" (gravastars) . They considered that independently of the matter that composes 
the gravitating object, as some critical limit is reached, the interior spacetime suffers a 
gravitational Bose-Einstein condensation [6]. The general relativity is not valid in the 
zone of coexisting phases, but it is valid macroscopically out of the transition zone [5], 
[6]. The net effect of the condensation is that the mean value of the vacuum stress- 
energy tensor changes from a nearly zero value to a non-zero value. The mean value 
of the vacuum stress-energy tensor has the form: < T^ u >=< p vac > g^, [9], [10J, 
which behaves like a cosmo logical constant term with A = 8nGc~ 2 < p vac >. This is 
equivalent to dark energy of the cosmological constant type, and behaves like a fluid 
with an equation of state: P = —pc 2 . The MM model is a static, spherical symmetric, 
five layer solution of the Einstein equations. It was built cutting and pasting three 
different exact solutions to the Einstein field equations. It has an interior de Sitter 
spacetime, while the exterior is Schwarzschild. The intermediate zone is a thick shell of 
matter that does not violate the strong energy conditions [UJ. The three solutions are 
matched by two thin anisotropic layers with distributions of surface tension and surface 
energy density. The solution does not possess an event horizon, but has a compactness 
that is close to one (the black hole compactness). 

Several researchers have analyzed the gravastar solutions using semi-analytic or 
analytic methods. For example, Cattoen et al [12] proposed a thick shell anisotropic 
gravastar model which has the advantage of smoothing out the large pressure jump 
present in the original MM model. Dymnikova [10] found new analytic solutions of the 
gravastar type, G-lumps, and A-black holes. Related new solutions of nonsingular black 
holes were also found by Mbonye and Kazanas [13J. Their solutions have an interior 
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de Sitter geometry containing matter. Chirenti and Rezzolla [2] found thick and thin 
shell gravastar solutions and studied their stability. Chan et al. [15] studied anisotropic 
phantom energy stars; while Lobo [16] studied gravastar solutions with a quintessence- 
like equation of state. Nicolini et al. [7], [8], obtained "mini-gr avast ar" and regular 
black hole solutions. 

Due to the large compactness of the gravastars it could be difficult to distinguish a 
gravastar from a black hole. Some arguments against the existence of gravastars, that 
could be verifiable through observations, were given by Broderick and Narayan |17j . On 
the other hand, positive detection features must include high energy particles emitted 
from the surface of the dark energy stars, as pointed out by Barbieri and Chapline [18J. 
The model of dark energy stars, or gravastars, is an interesting alternative to black holes 
and it is worth to explore its characteristics. 

A fluid of the cosmological constant type has necessarily a constant density in the 
absence of matter or other fields [10]. The reason is that for a source of the type 
Tpv = (87iG)~ 1 Ag^ u , the conservation equation T^ u - v = 0, implies a constant A. Thus, 
the dark energy must be coupled (somehow) to matter in order to trace the matter 
condensation. Another possibility is to consider a pure inhomogeneous anisotropic 
vacuum like term [6], [10]. In this work these two possibilities are combined, because 
the quantum vacuum phase transition must be driven by the matter during the stellar 
gravitational collapse. Thus, a set of compact objects made of fermionic matter coupled 
to inhomogeneous, anisotropic, dark energy is studied in this paper. 

The obtained solutions form a one-parameter family, where the parameter is the 
coupling constant between the dark energy and matter. The solutions converge on a 
model of the gravastar type as the parameter approaches a finite value. In addition, 
a new method is reported for the numeric integration of the anisotropic Tolman- 
Oppenheimer-Volkoff equations with cosmological constant (A-TOV) over the whole 
star. 

In the gravastar and dark energy models the shell is located above the event horizon. 
The exact position of the shell can not be given as a boundary condition of the A-TOV 
equations (without relaxing some other boundary condition), because the coordinates of 
the Schwarzschild radius are determined by the full solution. In this work, the position 
of the thin shell is obtained self-consistently as a result of the numeric integration. The 
mass-radius and structure of the stable models are estimated and compared with the 
Tolman-Oppenheimer-Volkoff neutron stars. 

The paper is organized as follows: in section [2] the Einstein equations and the 
notation are introduced. The A-TOV anisotropic equations are derived. In section [3] 
the equation of state for the matter and dark energy are given. In this section, the 
surface tension is obtained as a function of the anisotropy. The junction and boundary 
conditions for the A-TOV equations are discussed. In section HI the numerical algorithm 
is explained. In section the results are discussed. The paper ends in section [6] with 
some final remarks. 
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2. Einstein equations 



W u _ _ g „u R = , (1) 



The Einstein equations are: 

2 y c 4 

where i? denotes the scalar curvature, W u is the Ricci tensor, and is the energy 
momentum tensor. 

Assuming spherical symmetry the line element in standard coordinates [19] is: 

ds 2 = -e v{r) dt 2 + e x{r) dr 2 + r 2 {dd + sm 2 6d<j) 2 ) . (2) 

The energy-momentum tensor is composed of matter with mass-energy density 5 and 
pressure P, plus dark energy with density pa e , radial pressure Pd e (r), and tangential 
pressure P de (t)- 

T° = 6c 2 + Pde c 2 (3) 

Tl = ~ (P + Pde(r)) ( 4 ) 

Tl = -(P + P Mt) ) = Ti (5) 
I; = = (6) 

In terms of the (variable) cosmological constant the dark energy density is: p de = 
Ac 2 /87rG. The dark energy radial pressure is proportional to the dark energy density: 

Pde(r) = -Pde (? . (7) 

The variable dark energy density is assumed to be proportional to the mass density 
Pde — Qf p, where a is a non-negative constant (see section [3]). To sum up, the subindex 
"de" (dark energy) will not be used from now on: Pd e (r) = P r and Pd e (t) = Pt, except 
for the dark energy density pd e - 

The components of the Einstein equations are [20] . [T9] : 

T ° : 8vr(<5 + p de )/c 2 = (8) 



' X 1 \ 1 

7-h + H 



(9) 



Tl : -8n(P + P r )/c 4 = (10) 

An integral of equation (jSJ) is (see [20J): 

x 2mG , , 

e- X = l r , 12 

The mass-energy up to the radius r' is: 

m{r')=A7t f {5 + p de )r 2 dr, (13) 



o 
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In the special case of a constant dark energy density: m = m' + \{^-) r ^ , where m' is 
the integral (fT3l) performed over the mass-energy density of the matter alone. 
Subtracting Eq. ( flUl) from Eq. (jSJ) it is 

^^^ = -^(5 + P/c 2 ). (14) 
r ar c K ' 

At the surface of the star the right side of the equation is zero, so A + v is independent 
of r. To get an asymptotic flat solution it should be A, v — > (so A + v — > 0) as r — > oo, 
thus: 

A = — f for r > r s , (15) 

implies gr# = gr" 1 for r > r s , where r s is the surface radius. It can be seen from equation 
ffTiZj) that m(r) —>■ 0, as r — > 0, in order to get a regular metric at the center. So, an 
integration constant in equation (fTBj) was set to zero. 



2.1. Hydrostatic equilibrium structure equations 

The equations above give the metric coefficients as a function of the radial coordinate. 
The equation for the pressure profile can be obtained from the stress-energy tensor 
conservation. The Einstein equations and the contracted Bianchi identities give: 

r„?„ = o. (16) 

In hydrostatic equilibrium: 

= P, + ^(5c 2 + P) + -(Pr -Pt) = 0, 

' p a r 

where the notation a = e v , is used, and P = P + P r is the total radial 
pressure. Rearranging the terms the anisotropic A-Tolman-Oppenheimer-Volkov (A- 
TOV) equation is obtained: 

jp (mG | 4ttG p 3\ 

6/ -(5c 2 + P) ( 2 ■■' ' 



c?r r(r — 



2mG^ 



+ 2^, (17) 
r 

where AP = P t — P r , is the anisotropic term. This equation and its solutions were first 
studied by Bowers and Liang [2T] . Only for completeness, it is possible to rewrite the 
equation above for matter and uncoupled dark energy: 

alP AP 

-1~ = 2 + 

ar r 

(*!G + 4KGp 3 _ §gG 3) 

1 + j r(r - ^ - Prfe r3) ' 
or as function of the cosmological constant A = 8irG pd e / c 2 : 
alP „AP 
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This equation reduces to the well known Tolman-Oppenheimer-Volkov [22J equation 
when AP = 0, and A = 0. In this paper, equation (fTT|) is numerically integrated. 

3. Equation of State 

The equations obtained above must be supplemented with equations of state for the 
matter and for the dark energy. 

3.1. Equation of state for the matter 

A gas of neutrons at zero temperature has been considered here. This has led to a direct 
comparison of the results with the Oppenheimer-Volkoff neutron star model [22]. The 
energy per unit mass of the matter is given by |20j : 

4 5 1 

mtc 1 



6 c 2 



X 



Vl + x 2 (1 + 2a; 2 ) - ln(z + Vl + x 2 ) 



h 3 8tt 2 

This gives the total energy of the gas, including the rest mass energy density p, measured 
in [erg cm -3 ]. The Fermi x parameter is: x = p/m n c, where p = is the momentum 
of the particles, and n is the number density of neutrons. It can be written: 

,1/3 

x 



with: 



PqJ 



4 3 
7Th C 

P° = ^3 = 6 - 106 x 1015 fe cm ~ 3 ]- (2°) 



The pressure is given by [20] : 

' (21) 



4 5 1 

m„c 1 



h 3 8tt 2 

-2" 



x Vl+x 2 (2x 2 /3 - 1) + InO + Vl + x 2 ) 



measured in [dyn cm - ] . Note that usually the energy density is written as the rest mass 
density plus the internal energy density of the gas: 5 = p (1 + e/c 2 ) (see [20J). 

In this work the matter is assumed isotropic, while the dark energy is assumed 
anisotropic. 

3.2. Equation of State for the dark energy 

The dark energy density is assumed to be proportional to the rest mass density: 
Pde — oi p, where a is a non-negative proportionality constant (or coupling parameter). 
The dark energy suffers a phase transition at a certain critical mass density p c . This 
is another arbitrary parameter, but without loss of generality, a fixed critical density: 
p c = 2 x 10 14 g cm -3 has been chosen. If the mass density is higher than this critical 
value the dark energy switches to a non-zero value. Thus, the radial dark pressure is: 

P r ={-"^ T"-" C (22) 
tor p < p c 
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The tangential dark energy pressure and the radial dark energy pressure differ only 
at the shell of phase transition: P r ^ P t , at p = p c , and are determined by the field 
equations. 




o 5 10 15 

Radius [km] 



Figure 1. Pressure profile for each star. The pressure tends to a constant negative 
value inside the core as the critical parameter (a c = 0.5533) is approached. The 
pressure is measured in [gem -3 ] and normalized by the factor 10 13 c 2 . 



3.2.1. Anisotropy and surface tension There is a surface tension at the interface of 
phase transition. This is expected, in general, at the interface between two substances. 
In the absence of matter equations for pure anisotropic, inhomogeneous, dark energy 
[TO] are obtained from equations i fPTj) and ( 1221 ): 



P — — 2 

« = * + P3) 

A surface tension can be defined as: 

a = (AP) 8r, (24) 

measured in dyn/cm. Here 8r is the width of the shell, which is expected to be of 
Planckian thickness lp ~ 10~ 33 cm. With these definitions, equation (1231) can be written 
as: 

[P r ] = 2<7/r, (25) 

where [P r ] = P r + — P~, is the jump in the dark radial pressure, P^ is the radial pressure 
evaluated above/below the shell, respectively, and 2a /r is evaluated at the shell's radius 
r. Equation (!25l) is analogous to a Young-Laplace equation for the spherical interface 
between two substances [21]. The anisotropy (surface tension) compensates the large 
dark energy pressure gradient. From Eq. (1231) above, it follows that a positive surface 
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tension implies a larger radial dark pressure from above the thin shell yj. There is a 
net radial pressure compressing the shell. In this case, the tangential dark pressure 
is larger than the radial one on the shell (see Eq. 1241) . On the contrary, a negative 
surface tension, gives a net radial dark pressure directed outwardly, and a tangential 
dark pressure lower than the radial one. 

The surface tension is not directly calculated (and it is not necessary) in the 
algorithm, but it could be found with an indirect calculation. It can be observed that 
it remains finite. This is inferred a posteriori from the numeric results, using the Israel- 
Lanczos [25] thin shell junction conditions. As observed by Mazur and Mottola [6], 
the coefficient g u of the metric is continuous, but its first derivative and the metric 
coefficient g rr are in general discontinuous, due to the discontinuity in the equation of 
state. The outwardly directed normal vector to the interface is n = ^fg~^-d r ^ and the 
extrinsic curvature is defined as K\ = V a n 6 , where a, b, are indices of the hypersurface 
coordinate system. The discontinuities in the extrinsic curvature as a function of the 
surface stress energy density r\ and the surface tension a are [6]: 



The first condition given by Eq. ([26]) is related to Eq. ([25]) . An assumed jump in the 
radial pressure gives a surface tension through Eq. (|23|) . while Eq. ( 12"6"1) gives a jump 
in the first derivative of g a (there is no jump in g rr , since it is assumed that rj — 0). 
Conversely a jump in the first derivative of gu is proportional to a surface tension, and 
this could be calculated after obtaining the numeric solution (see the next section). 
Equation (|25l) is valid at the shell even when the isotropic fermion matter is included. 
This is described in the next section. 

A word of caution is worth mentioning. It is possible to derive an expression, from 
equations ( |23i) and ( 122"1) . for the tangential dark energy pressure at the shell, analogous 
to that obtained by Israel and Poisson [I]: 



where a is the radius of the shell. They pointed out that this expression becomes 
singular at the event horizon. In fact, as they showed, a change in variables gives 
5(r — a) = [<?rr(a)] 1 ' /2 3( s ), m terms of the proper distance s from the shell. So, if r = a 
is light-like the pressure becomes singular, "even when considered as a distribution" [4]. 
This means that the thin shell of phase transition cannot be located exactly at the event 
horizon, but can be located elsewhere. Although this conclusion was originally obtained 
by a precursor idea of the dark energy stellar models, it is of course, applicable to the 
models studied here. 

| In fluid dynamics [24] the Eq. ([25]) is commonly written — [P r ] = 2a' /r , with a' = —a. Thus, what 
is called a positive surface tension <j here is a negative surface tension in the textbooks, and vice versa. 




(26) 




(27) 
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Radius [km] 

Figure 2. Mass density profile for each dark energy star. The mass density is 
not smooth at the phase transition shell. The density is measured in [gem -3 ] and 
normalized by the factor 10 13 . 



4. Numeric algorithm 

The discontinuity in the dark pressure is managed as follows: the radial derivative of 
the total pressure in Eq. (|T7|) is split in two terms, 

dP , (5P r 8P\ . , 

— = lim — ± + — , 28 
dr Sr-^o \ 5r 5r J 

taking the limit with 5r tending to zero. The difficulty to integrate the A-TOV equations 
over the whole star is that the first term on the right side of Eq. (12 8j) diverges at the 
shell: lim^^o ^\ r=a —> oo. This term is associated with the second term on the right 
side of Eq. f l27|) . The solution, proposed here, is to cancel that term with an appropriate 
amount of anisotropy on the right side of Eq. (fT7|) . As explained, this is equivalent to 
set a finite surface tension at the shell (see Eq. [241) . Now the set of equations to be 
integrated are equations (TT2"j) . (ITB1) . and (JT7J) (without the discontinuity), complemented 
with the equations of state described above. The boundary condition for the Eq. ffTTj) 
is P = at the surface of the star, and the initial condition, a prescribed mass density 
at the center. 

Briefly, the algorithm works as follows: first, it integrates the equation 

jp (jnG , 4ttG p 3\ 

f = + ' (29) 

from the center of the star, with a given central mass density. Next, it updates the 
values of P, P, p, 5, pd e , and m. If the mass density is above zero, it goes a radial step 
further and repeat all. If it is zero the code ends. Of course, in the practice 5r is never 
zero, but equal to the radial step of the integrator with 5r a. It has been checked 
that the results are robust respect to variations in Sr. 
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Note that the Eq. (I29p is different from the Eq. (1171) . The term dP r /dr was 
canceled by a surface tension a so that 5P r = 2a /r. 

Equations ( [121) . ( TT3T) . and ( fI7|) are integrated with a Runge-Kutta algorithm of the 
fourth order. The employed code HE05v2 is an update of the previous code HE05vl 
(see Ref. [20] for more details). 




4-1- Caveats 

Although the dark energy pressure derivative was canceled, it is not granted that the 
solutions will be smooth. The solving procedure described above eliminates the delta 
Dirac distribution from equation (fTTl) . but smoothness was not imposed anywhere in 
the remaining equations. As it will be shown, there are solutions which are not smooth 
at the transition shell. It is possible to integrate the non-smooth solutions performing 
the calculation at zero order precision on the shell, and at fourth order in the rest of 
the star. 

Another caveat is that the denominator of the Eq. (TT7T) is zero at the event 
horizon, and the equation becomes singular. A larger numeric resolution is needed 
as the compactness approaches one. This problem is inherent to the field equations, 
and of course, it is independent of the position of the transition shell. However, it also 
includes the issue pointed out above: the transition shell cannot be located at the event 
horizon. 

5. Discussions 

A family of solutions is obtained, depending on a, for each stellar central density. 
In order to analyze the structure of the stars, a model with central density: 
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Figure 4. The compactness function: r sc h/r, is plotted for several dark energy stars. 
As the parameter a is increased, the compactness at the surface tends to one, and the 
model approaches a regular black hole. 

p = 2.86 x 10 14 gem -3 , is singled out. The results of the calculations are shown in 
Figures [THU Figures [Q [2], and [31 show the stellar structure of the set of dark energy- 
stars. These stars are in the stable branch of the solutions (see below). The different 
solutions were obtained varying the coupling parameter from a = to a = 0.5533. 
The summary of possible solutions is given in Table (H presented together with some 
well known solutions. The type (f ) solutions shown in the Table [H are obtained setting 
a = 0, and correspond to the typical Tolman-Oppenheimer-Volkoff neutron star model. 
Black hole solutions are obtained with a > a c = 0.5533. The critical value a c is obtained 
by interpolation and cannot be calculated exactly [§[ as was explained above (see Sec. 
13.2. II and l4l). The black hole solutions are characterized by a surface radius smaller than 
or equal to the Schwarzschild radius, or a larger or equal to one compactness. These are 
indicated as type (a) and (b) solutions in Table [H The critical a c parameter depends 
on the particular family of solutions considered, i.e.: on the central mass density. 

Figure [1] shows the pressure profile for each dark energy star of the family. The 
pressure tends to a negative constant inside the dark energy star as the critical parameter 
is approached. The pressure is positive and decreasing with radius outside the thin 
shell. Figure [2] shows the mass density for each dark energy star. The density tends 
to a constant value, inside the shell, as the critical parameter is approached (a — > a c ). 
The density profiles are non smooth for this family, because there is a surface tension 
at the thin shell of phase transition (see Eq. 1261) . This is expected in general, since the 
matter density need not be smooth across the shell. 

§ In fact the interpolated value is a c = 0.553317, with a precision of the order of 10 cm in the radial 
step. 
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Table 1. Summary of solutions. 



type w a 



matter EOS layers Anisotropy converges to 



a < 1 > a c (a < a' c ) neutrons 

b > 1 > a c (a > a' c ) neutrons 

c 1 = (a' c < a c ) neutrons 

d < 1 < min(a' c , a c ) neutrons 

e > 1 a' c < a < a c neutrons 

f - neutrons 



3 
3 
3 
3 
3 


3 
5 



Yes 
Yes 
Yes 
Yes 
Yes 
No 
Yes 
No 
Yes 



Black Hole 
Black Hole 



TOV star [2J 

anisotropic star, Bowers-Liang [21 
DE star, Chapline et al. [5] 
Gravastar, Mazur-Mottola [6 



Gravastar, cosmological const, as r — > 



Gravastar 
Gravastar 



P = P /3 



1 - P = const. 

1 - P = const 



In general, the total pressure (effective equation of state) can be written as: 
P = —w pc 2 , 



(30) 



with 



w 




) 



(31) 



The value of w is bounded from above by a, and it depends on (a, p). For the analyzed 
model w = 0.5287 for a = a c . Thus, the total pressure is P > — pc 2 , for the considered 
set of dark energy stars. Every family of solutions contains at least a black hole, i.e.: 
the critical parameter a c is finite for all families of solutions. This was obtained by 
inspection of the numerical solutions. 

For a given value of the central density (for a given family of solutions), there exists 
an a = a' c , so that the effective equation of state is of the cosmological constant type 
at the center of the star, i.e.: P = —pc 2 (do not confuse a' c with a c ). If a c > a' c the 
sequence of solutions, as function of a, converge on a dark energy star of this kind (type 
c) before converging on a black hole. On the contrary, if a c < a' c , the sequence converges 
first on a black hole. 

The critical parameter a' c can be obtained exactly from the neutron gas equation 
of state. For example: the Fermi parameter is x = 

(p/Po) 1/3 = 0.3604 for the solutions 
considered here (with central density p = 2.86 x 10 14 gcm" 3 ). With w = 1, using Eqs. 
flM]) and (TJU), is obtained: 



that gives a' c = 1.0248. The same calculation can be done for any other central 
density. The value of a c must be obtained numerically, i.e.: it cannot be obtained 
without knowledge of the complete sequence of solutions. In the present case a c < a! c 
so, the solutions converge on a black hole before reaching a solution with an effective 
cosmological constant equation of state, at the center of the star. The parameter space 
was swept looking for another family of solutions that contains a type (c) gravastar. 
This is not a trivial task, because w must be calculated for a = a c , while a c vary 



a ' c = 1 + - \x Vl + x 2 {2x 2 /3 - 1) + ln(x + Vl + x 2 )} /x 3 , 



(32) 
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between the different families. So, a complete set of simulations had to be performed 
for each central density. It was found that w(a c ) increases with the central density 
of the models. In particular, a gravastar model of type (c) was found with a central 
density of p = 4.92 x 10 16 gem -3 . In this case, the critical parameter for the type 
(c) gravastar is a' c = 1.416, while the critical parameter for a black hole solution is 
a c = 1.5192. The compactness of the gravastar is 0.9079. The structure of this star is 
shown in Fig. [51 Unfortunately, this model is in the unstable branch of the solutions. 
Of course, as explained above, it is possible to find solutions -stable or unstable- with 
larger compactness, but with w 1. 

There is another type of solutions which are closer to the MM gravastar model: the 
solutions with a — > a c , which have a density profile nearly constant in a large portion 
of the stellar interior (see Figs. [Hand [2]). The models with a — > a c -but not reaching 
a c - have a large compactness and are not contained in a black hole. In order to have 
an inner de Sitter metric, the matter density must be exactly zero at the core; this can 
be checked from Eqs. (1121) . (TT3T) . and (114j) . However, the numerical limit of the models 
as p, p c — > do not converge to the MM model: the density and pressure profiles keep 
their shape roughly, but the radius of the dark energy stars increase continuously as the 
density is lowered. It means that the radius of the stars tends to infinity as the central 
density tends to zero. The limit spacetime in this case is Minkowski (vacuum). But 
there is another way around to obtain an MM gravastar, with an inner de Sitter metric, 
as the limit of the numerical model: it is to change the order of the limits. First kept 
fixed pde =constant in the obtained solution for a — > a c . This is equivalent to leaving 
the dark energy and shell radius frozen. After that, take the limit p — > 0, to obtain an 
inner de Sitter spacetime. In the process the inner mass will be overestimated, so the 
result is a perturbed MM gravastar. These solutions were not included in the Table HJ 
but can be derived from type (d) or (e) solutions. It is not possible to know how this 
gravastar will evolve with the present methods, so this issue will not be further discussed 
here. 

The compactness is defined as the quotient between the Schwarszchild radius and 
the radius of a sphere. It is a function of the radius and the parameter a. Figure 0] shows 
the compactness function for the considered family. It can be seen that the compactness 
of the star surface tends to one as a — > a c . If a > a c the compactness is larger or equal 
to one, and the shell is contained within a black hole. The solutions with a > a c are 
black holes, indicated as type (a) or (b) in Table [U 

Figure [6] shows the mass-radius relations for dark energy stars, with a coupling 
parameter that varies from a = to a = 0.53. There are two maxima in the mass- 
radius curve. The first one is at the density of phase transition, where the curve is non 
smooth. The second maxima is due to the matter equation of state, like in the usual 
TOV neutron stars. There is a valley in the curves because at low masses the repulsive 
effect of the dark energy is more important than its weight. However, at larger masses 
the weight overcomes the negative pressure effect and the masses of the dark energy 
stars become larger than fiducial TOV stars with the same radius. 
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Detail of the pressure and density profile with a 
gem" 



a' and with central 



density p — 4.92 x 10 16 g cm 3 . In this case the magnitude of the total radial pressure 



at the center is: P — —pc 2 . The physical units of the density and pressure (xc 2 ) are: 
[le!6] gr cm -3 , the scale is linear. 



The dark energy star of the family is initially (at a = 0) on the stable region near 
the first maximum, but its position on the mass-radius curve is displaced to the left, to 
lower radius, as a increases. 

The criteria dM/dp > 1 (dM/dp < 1) is a necessary condition for hydrodynamic 
stability (instability), respectively, but it is not a sufficient condition |23j. However, it is 
possible to apply a Wheeler's M(R) analysis to the mass-radius relation. The method 
was devised by Wheeler, and the precise conditions under which it is valid, were given 
by Thorne (see Bardeen et al. [26]). It is useful, in the present case, to obtain the 
number of unstable radial modes for each dark energy star in hydrostatic equilibrium 
without calculating the eigenfunctions and eigenfrequencies of the radial perturbations. 
The condition for the applicability of the method is satisfied: the microscopic equation 
of state represents cold catalyzed matter, i.e.: no changes in nuclear composition are 
considered [JJ]. The method is applied as described in Ref. [26]: the graph of the 
mass versus radius, is parameterized by the central density of the configuration (for a 
particular a value). As we move in the direction of increasing density all configurations 
remain stable until the first maximum value of M is reached. At the first extremal point, 
the fundamental mode of radial oscillation becomes unstable. At the second extremal 
point, another mode of radial oscillation changes stability. The direction of stability 
change depends on the shape of the M(R) curve at that extremal point. If the curve 
bends clockwise, with increasing density, then one previous unstable mode becomes 
stable. On the contrary, if the curve bends counterclockwise, then one previous stable 
mode becomes unstable. The same criterion is applied for each succeeding extremal 

|| For a more realistic equation of state the Wheeler M(R) analysis is applicable if the adiabatic index 
satisfies a certain criteria 26 . 
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It is generally assumed that a cold star can be dynamically unstable against 
nonradial perturbations only if its lowest radial mode is also unstable (see [27]). Thus the 
radial stability analysis reveals the absolute stability and instability of the dark energy 
stars. According to this analysis, the stability and instability regions for a = 0.53 are 
indicated in the Fig. [6] with the letters "s" and "u" respectively. There appears an 
unstable band between two stable ones, starting at the density of phase transition and 
ending at the bottom of the valley. This is analogous to the unstable band between the 
white dwarfs and the neutron stars (see [26]). 




Radius 

Figure 6. Mass radius relation for a set of solutions. The curve with a = corresponds 
to normal TOV neutron stars. 



6. Final Remarks 

In this work, the A-TOV anisotropic equations were solved for compact objects 
composed by a Fermi gas coupled to dark energy. The dark energy traces the matter 
distribution at high energy density, with a strength depending on the parameter a. 
The main difference is that, in the original dark energy stellar model, the core of the 
star is de Sitter and the shell is located close, or replacing, the event horizon. In the 
solutions found in this paper, the position of the shell is calculated from the model and 
not imposed as a boundary condition. Moreover, unlike other gravastar models, the DE 
stars considered here contain matter in its core (an exception is the model of Ref. |13j). 

A simple numerical method was devised to integrate the stellar structure equations 
from the center of the star to its surface, including the shell. The compactness of the 
stars, their structure and mass-radius relations were calculated for comparison with the 
well known theoretical TOV model. 

% This description of the Wheeler analysis is based on the paper of Bardeen et al. [26] . Another clear 
exposition is given in the reference [27] • 



Ghezzi C. R. 



16 



It is found that the dark energy plus anisotropic stress helps to stabilize a neutron 
star, but there is yet a maximum mass for stable DE stars. 

It has been shown that the effective equation of state is: P = —wpc 2 , with w = 1, 
w < 1, or w > 1, depending on the density and the assumed value of a. There is 
a critical parameter a' c so that the total pressure, at the center of the star, is of the 
cosmological constant type P = —pc 2 . 

There exists, also, a critical parameter a c which separates dark energy stars from 
regular black hole solutions. 

The mass-radius curves for the dark energy stars have a maximum mass, signaling 
the unavoidability of black hole formation for large enough mass. It is unlikely that 
some other EOS for the matter will avoid the BH formation. The only possibility left is 
a different interaction law between dark energy and matter, for which we have no clues 
at present. This issue must be further studied. 

One of the characteristic features of the mass-radius relation for the dark energy 
stars is that there appears an instability gap between the two stable bands, which do not 
exist for TOV neutron stars. It seems probable that these unstable stars migrate to the 
stable branch at larger radius. The detection of this gap could imply a large interaction 
between dark energy and matter, and could be a clear signature of the validity of the 
model. However, its negative detection does not rule out the existence of dark energy 
stars. In this case, their identification through the mass-radius curve could be more 
difficult. 
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